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Abstract 


In  this  paper  we  discuss^  the  regularity  of  the  eigenfunctions  of  eigen¬ 
value  problems  with  piecewise  analytic  data  and  the  approximation  of  the 
eigenvalues  and  eigenfunctions  of  such  problems.  A  detailed  and  systematic 
numerical  study  of  these  approximations  is  presented,  together  with  an 
analysis  of  the  numerical  results  in  light  of  the  theoretical  results.  The 
specific  aim  is  to  assess  the  reliability  of  the  theoretical  results  -  which 
are  of  an  asymptotic  nature  -  as  a  guide  to  practical  computations  -  which 
may  take  place  in  the  preasymptot ic  phase  -  and  to  look  for  characteristic 
features  of  the  numerical  results  which  are  not  completely  explained  by  known 


theoretical  results. 


Regularity  and  Numerical  Solution  of  Eigenvalue  Problems 
with  Piecewise  Analytic  Data 

1.  Introduction 

The  purpose  of  this  paper  is  to  discuss  the  regularity  of  the  eigenfunc¬ 
tions  of  eigenvalue  problems  with  piecewise  analytic  data  and  the  approxima¬ 
tion  of  the  eigenpairs  of  such  eigenvalue  problems  by  the  finite  element 
method,  and  to  present  the  results  of  a  detailed  and  systematic  numerical 
study  of  these  approximations. 

[1]  contains  an  analysis  of  the  regularity  properties  of  the  solutions  of 

elliptic  boundary  value  problems  with  piecewise  analytic  data.  Specifically, 

it  is  shown  that  the  solution  of  such  a  problem  belongs  to  a  countably  normed 
2 

space  B„(n).  In  [2]  this  fact  is  used  to  show  that  the  h-p  version  of  the 

p 

finite  element  method  has  an  exponential  rate  of  convergence.  (For  a  survey 
of  the  basic  results  on  the  h-p  version  of  the  finite  element  method  and  for 
relevant  references  we  refer  to  [3]  and  [4]. )  In  this  paper  we  make  a 
parallel  study  of  the  regularity  and  approx i mabi 1 ity  of  the  eigenfunctions  of 
eigenvalue  problems  with  piecewise  analytic  data.  We  then  discuss  the  impli¬ 
cations  of  this  approximabi 1 ity  for  the  approximation  of  eigenvalues  and 
eigenfunctions  by  the  finite  element  method. 

Section  2  contains  background  information,  Section  3  describes  the  eigen¬ 
value  problems  we  treat,  Section  4  presents  a  regularity  result  for  eigenfunc¬ 
tions  in  terms  of  countably  normed  spaces,  Section  5  surveys  (abstract) 
results  on  eigenvalue  and  eigenfunction  approximation  by  finite  element 
methods,  and  Section  6  reports  the  results  of  a  detailed  and  systematic 
numerical  study  of  eigenvalue  approximation  by  finite  element  methods  and 
presents  the  conclusions  of  this  study. 


2.  Notation  and  Preliminaries 
2 

Let  Q  c  IR  be  a  simply  connected,  bounded  domain  with  boundary  r  =  6Q, 
M 

and  assume  that  T  =  UP.,  where  P.  ,  for  i  =  1 . M,  is  an  analytic 

i  =  l  1  1 

simple  arc  connecting  the  vertices  A.  .  and  A.  (A^  =  Aw).  r,  will  denote 

i-l  i  0  M  i 

T.  -  {A  .,A.}.  The  interior  angle  at  A.  will  be  denoted  by  o>.  .  This 

li-li  l 

notational  scheme  is  shown  in  Figure  2.1. 


A3 


Figure  2.1.  Notational  Scheme  for  n  and  F  =  5fi. 

We  will  assume  0  <  on  <  2tt.  =  2rr  corresponds  to  a  slit  in  Q.  If  w.  < 

2n  for  all  i,  we  will  call  Q  a  Lipschitz  domain  (T  will  be  a  Lipschitz 
curve  in  this  case).  If  the  arcs  r  are  straight  lines,  we  say  Q  is  a 
straight  polygon  or  simply  a  polygon;  otherwise  we  wili  refer  to  Q  as  a 
curvilinear  polygon.  We  let 

r°  =  UP., 

leV  1 


where  T>  is  a  subset  of  {1 


M} ,  and  let 


boundary.  We  will  often  consider  the  special  domains 


S  =  S(w,5)  =  {(r.Q)  :  0  <  r  <  5,  0  <  0  <  u} 


where  (r,0)  denotes  the  polar  coordinates  of  x.  In  connection  with  S  we 
introduce  the  notation 


r^w.5)  =  { ( r,  0 )  :  0  <  r  <  5,  0  =  0}, 


5)  =  {(r,0)  :  0  <  r  <  5,  0  =  u) . 


By  H  (Q),  for  m  >  0  an  integer,  we  denote  the  standard  Sobolev  space 
of  functions  u  for  which 


(2.  1) 


Hm(n) 


=  £  IID“ul 


L2(t2) 


a  d  I  “I  u 

where  D  u  =  - ,  with  a  =  (a  a  ),  a ,  a  >  0  integers,  and 

Sx^Sxf2  12  12 

a  +  a  .  We  note  that  ||u||  =  l|u||  ,  ..  Let 

H  (Q)  L2U  ; 


H^(n)  =  {u  :  us  u  =  0  on  f°} ; 


on  H  (Q)  we  can  use  either  |iu|| 

H1  (0) 


(  I 


r.(x)  =  dist  C  x,  A  . )  =  I  x  —  A  .  |  , 
J  J  J 

r.(x)  =  mi n( 1 , r  .  ( x) ) , 

J  J 


r.  <•.  -r.  V.  v'.-'  •  -\  w-  -V--.*'- ■r. 

v  V  '/  V  V  /  V  *  ■*.-  '  j  v  *  *  ■  m  * 


(3  =  . /3m),  where  0  <  (3.  <  1, 


k  =  an  integer, 


and  define 

M  *  8  +k 

4>  (x)  =  IT  (r.  (x)  r  1 

P+k  i=l  1 

If  n  =  S,  we  let  $  =  r^.  By  H^’^fl),  with  m  and  £  integers  satis- 

P  H 

fying  0  <  £  <  m,  we  denote  the  weighted  Sobolev  space  characterized  by  the 


'LL 


(2.2)  ||uirm 

H^(n) 


I  z 

k=0  |a|=k 


I*  lAii'o  •  e  =  °’ 

P  K  H  (0) 


H„  (Q)  will  also  be  denoted  by  L0(fi). 
p  P 

We  will  also  use  the  countably  normed  spaces 


S@‘°)  -  <“  :  “  6  “e,!(n>-  lVw“‘,lso(11I 


\r  —  P 

<  CdK  (k-m 


(2.3) 


for  | a |  =  k,  k  =  £,£+!,...,  where  C  >  1  and 


d  >  1  are  independent  of  k> 


■2(n)  =  {u  :  u  e  H?’2(n), |Dau(x) |  <  Cdkk! (*  1 (x) )' 


(2.4) 


for  I  a |  =  k,  k  =  1,2 .  where  C  >  1  and 


d  >  1  are  independent  of  k} . 


In  the  case  when  Q  =  S(cj,5)  we  let 


4 


3rai58a2  rai6a2 


jn  t 

and  then  define  the  spaces  (S),  with  0  <  l  <  m,  and  0  <  (3  <  1,  which 


are  characterized  by  the  norm 


(2.5)  ||u!! 


,2 


(lull  „  + 

"  "  p-\ 

H  (S) 


£  £  iir^1  v'Hrurft  .if  *  >  i 


^ai-^+p.a  2 


k=£  !  a  j  =k 


H°(S) 


I  I 


k=0  |a|=k 


|rai+02)au!|2  ,  l  =  0. 

n  (S) 


,0,0, 


K  ’  (S)  will  also  be  denoted  by  if  ( S ) .  Note  that  L0(S)  =  if _ ( S ) .  Further- 

p  p  p  p 


.1 


more,  we  let  H  (S)  be  the  space  characterized  by  the  norm 


Hull2.  =  Hull2  +  V  ||rai  ^“ull2 

X (S)  HU(S)  HU(S) 

|a|=l 


Note  that  this  is  just  [|u||  expressed  in  polar  coordinates. 

H1  (S) 


We  now  state  some  lemmas  we  will  use  in  the  sequel. 


Lemma  2.1.  The  spaces  H^’J(S)  and  are  equal  for  j  =  0,1,2, 


i .  e.  ,  the  norms  ||  u||  . 


J(S) 


and 


J(S) 


are  equivalent  for  j  =  0,1,2. 


(3 


This  result  is  obvious  if  j  =  0,1.  For  a  proof  in  the  case  j  =  2,  see 
Lemma  1.1  of  [ 1 ] . 


Lemma  2.2.  Suppose  0  <  l  <  2  and  k  >  i.  Then 


(2.6) 


H*W_pDaull  n  ^  Cdk~*(k-£)!  for  any  a 
1  H  (S) 


with  |a| 


1 


if  and  only  if 


(2.7) 


'  ?  _2ao 

|Da  u|r  ^<t>)3+k,_£rdrde 


<  Cdk'  £( k' -£) !  for  all  a' 


with  l  <  |a' t  =  k'  <  k. 


For  the  proof  of  this  result  see  Theorem  1.1  of  [1  ], 


Lemma  2.3.  Let  u  e  H  (S(u,<5))  satisfy 


(2.8) 


-Au  =  f  on  S(u,5) 


and  suppose 


(2.8a) 


u  =  0  on  u  1^2- 


(2.8b) 


u  =  0  on  f, 


~  =  0  on  r 
6n  2 


(2.8c) 


o  u  pi  pi 

—  =  o  on  r  u  r  . 

on  12 


(2.9)  ||u|!  <  C  ||f ||  -n  +  Hull 

VT'  (S(w,5/2))  »■  (3(  (  U  H  '  1(S(u,5)-S(w,5/2)  ) 

p  P 


for  any  [3  satisfying  0  <  (3  <  1  and,  in  addition,  1  -  n/o)  <  ft  in  cases 
(a)  and  (c)  and  1  -  rr/2w  <3  in  case  (b). 


For  the  proof  of  (2.9),  see  the  proof  of  Theorem  2.1,  specifically  the 


proof  of  inequality  (2.44),  in  [1], 


Lemma  2.4.  2L(n)  c  c  20  (fi),  for  any  e  >  0. 

--  p  p  p+c 


For  the  proof  of  this  result  see  Theorems  2.2  and  2.3  of  [5 


3.  The  Eigenvalue  Problem 


Let  fi  be  a  straight  or  curved  polygon,  and  consider  the  selfadjoint 
eigenvalue  problem 


(3.  1) 


(Lu)(x)  =  A(Mu)(x),  x  =  (x^.x^)  6  Q 
(Bu)(x)  =  0,  x  e  T, 


where 


(3.2) 


(Lu)(x)  = 


Z  a^hj(x)  &<x)]  *c,x)u’ 


i ,  j=l 


with  a.  .(x)  =  a..(x)  and  c(x)  >  0  analytic  functions  on  Q  satisfying 
ij  Ji 

(3.3)  ^a.J(x)Ci?J  >  ?1(^+C2]  V  x  e  n  and  V  (f;^),  with  ^  >  0, 

i.  j 


where 


(3.4) 


(Mu)(x)  =  d ( x ) u ( x ) , 


with  d(x)  an  analytic  function  on  n  satisfying 


(3.5) 


d(x)  >  >  o. 


and  where 


u(x),  X  e  T 


(3.6) 


( Bu) ( x)  = 


3u,  ,  V  3u  ,  s  _1 

^ — ( x )  =  -  >  a.  .n.  -s — (x),  x  €  T  , 
L  IJ  J  3x  = 

i.J 


with  n  =  n(x)  =  (n^.n^)  denoting  the  exterior  unit  normal  to  V  at  x. 
we  define  the  bilinear  forms 


(3.7) 


a(u, v)  = 


.  ,3u  3v  "I  . 

a.  . ( x ) ^ —  ^ —  +  cuv  dx 
l.i  ax.  3x  . 

u  1  j  J 


(3.8) 


b(u,v)  =  duvdx 


on  Hq(H),  then  (3.1)  has  the  weak,  or  variational,  formulation 


(3.9) 


Seek  A  e  R,  0  *  u  e  HQ(Q)  satisfying 
a(u,v)  =  Ab(u,v)  V  v  e  Hq(Q). 


(3.9)  has  eigenvalues 


0  <  x:  <  x2  <  .  .  .  /’+« 


and  corresponding  eigenfunctions 


U1,U2* 


satisfying 


a'W  » 


In  Section  4  we  prove  the  following  theorem  on  the  regularity  of  the 


eigenfunctions  of  (3.9). 


Theorem  3.1.  Suppose  Q  is  a  straight  or  Lipschitz  curved  polygon  (as 
described  in  Section  2).  Then  for  each  j  =  1,2,..., 


u.  e  2T(n), 
J  P 


where  p  =  (/3, , .  .  .  ,  (3..) ,  with  8.  satisfying 
1  M  1 


0  <  ft.  <  13.  <  1, 

l  l 

where  /3 .  depends  on  the  values  a  (A.)  and  the  angle  w.  . 

l  mn  l  l 

Remark  3.1.  We  have  assumed  Q  is  simply  connected.  Theorem  3.1  can  easily 


be  extended  to  cover  multiply  connected  domains. 


9 


Remark  3.2.  In  Theorem  3.1  we  have  assumed  that  a.  .,c,  and  d  are  analytic 

^  J 

on  Q.  This  condition  can  be  weakened  to  the  requirement  that  these  functions 

M 

are  analytic  on  Q-  U  A  ,  provided  we  impose  certain  growth  conditions  on 

i  =  l  1 

a.  ,.c.d.  and  their  derivatives. 


jW^WAiVAWWWWiV« «^^.^'^VTVrV^’TVTJTO’Ty»VrJTO"^>'’jr'jB--jr«jrr^r'jr'>-vr'jrrjrVTJrVf- 


4.  Regularity  of  the  Eigenfunctions 

In  this  section  we  will  prove  Theorem  3.1.  We  provide  the  complete 
details  only  for  the  case  when  Lu  =  -Au  and  fi  is  a  straight  polygon  with 
r°  =  T  (so  that  the  boundary  conditions  are  of  Dirichlet  type).  The  proof  in 
the  general  case  will  only  be  outlined  since  the  arguments  are  similar  to 
those  used  in  [1]. 

M 

Lemma  4.1.  The  eigenfunctions  of  (3.9)  are  analytic  in  Q  =  Q  -  U  A. . 

i  =  l  1 

For  the  proof  of  this  result  see  [6,  Theorems  5.7.1  and  5.7.1']. 

Proof  of  Theorem  3.1.  Assume  Q  is  a  straight  polygon,  and  suppose,  without 
loss  of  generality,  that  A^  is  located  at  the  origin  and  lies  along 

the  x-axis.  Choose  0  <  5.  <  1  so  that  S_  =  S(w.  ,5.)  c  £1  and  S_  n 

1  1  Oi  1  1  0[ 


I  FJ  =  a' 


See  Figure  4.  1. 


j*i. i+1 


Figure  4.1.  The  Domain  S„  . 


Consider  an  eigenpalr  (A,u)  of  (3.9)  and  suppose  j|u||  =  1. 

IT(n) 

Localizing  to  the  vertex  A^  we  let 


(4.  1) 


v.  =  r  u_k,  k  =  0,1,2,. 


Then,  using  the  identity 


k  k-1  2  k-2  2 

(r  u  k)  =  r  (r  Au)  k  =  -Ar  ir  u)  k,  k  >  0, 
r  r  r 


we  have 


(4.2) 


■Av  =  Ark  *(r2u)  k  in  S_ 
k  r  Si 


vk  =  0  on  r,u  r2. 


By  Lemma  4.1,  u  is  analytic  in  fl.  Thus  there  exist  constants  and  d^ 

such  that 


(4.3) 


'k«  1  £C<Jk.,k-0,l . 

*  <  W?’ 


Let  us  now  prove  that  there  are  numbers  C  and  d  such  that 


(4.4) 


I v  ||  <  Cd  k!  ,  k  =  0,  1,2 . 

VC'  (S  ) 

(3  6l/2J 


for  any  0  <  0  <  1  satisfying  1  -  n/u>^  =  f3^  <  f3.  We  will  show,  in  fact, 
that  (4.4)  holds  provided  C  =  C^Cq+X+I),  d  satisfies 


3C  A  C  A 

-7T-+  *  1/2, 

d  d2 


cociu+1)do  2CociA  C0C1A  c 

d  d  +  2  ~  2’ 

d 


where  and  are  the  constants  in  (4.3)  and  (2.9),  respectively,  and 

(3.  <  (3.  First  consider  k  =  0.  u  =  v  e  H^S.  ) ,  u  =  v_  =  0  on  f.  u  f„, 

1  U  Oj  U  1  Z 

-Au  =  Au  on  S  (cf.  (4.3)),  and  ||u||  =  1 ,  so  we  can  apply  Lemma  2.3 

5‘  H^n) 


12 


to  obtain 


°V.2f<;  , 


^•2'S6,/2» 


SCi(: 


A|Iu||  +  l!u|| 

*  *  ** (s«1-sa1/2) 


<  C,  max (1,57) 
1  l 


A||u||  +  ||u|| 

HlS  )  n  (S_  -S_  ._) 

Oi  0\  O^/d 


<  C  (A+l )  ||u|| 

Hx(n) 


<  C^(A+1),  provided  0  <  /3  <  1  satisfies 


1  -  n/u>l  =  /3.  <  0; 


thus  (4.4),  with  k  =  0,  holds  for  this  choice  of  C  and  d. 

Suppose  (4.4)  holds  for  k  =  0,1  ,...,£- 1,  where  £  >  1.  We  wi  1 1  shew  i t 


holds  for  k  =  £  and  conclude  by  induction  that  it  holds  for  all  k.  We 


begin  by  noting  that  v  e  H  (S„  ).  To  see  this  it  is  sufficient  to  show  that 

t  di 


I v  ||  <  oo.  We  easily  see  that 

H  (S61/2) 


'J  ,  =  II  ( l-£)v  +  r(v  )  || 

H!(S  )  g~1  g~1  r  Hl(S.  yo) 

0[/d  O^/d 


s  c(»Vi«Hi 

£  ch-i" , 


+  lr(vM)rl  i 


+  llr(v»  1  )_ll  n 

H‘(SS,/2>  H  (S5,/2> 


‘»r,Vl'p'r"rf.(s 


S-,/21 


-1. 


*,r  lr(w(-lV'e"„0„ 


i 


wvi*  vwvry'f  wvvi'w  o v*  v»  >r«  jr»jr»j^»>TjrMrpjf  ">'"3rr^r"#T>  v  v  ’^’V  -w.  ■V'.'W’V  'aw.vwv.v 


s  c[,v<-illHi 


+  Br(Vl)rl  0 

H  (S5i/2} 


*,r(V'->l-"„o,s 


5 1  /2 


+  H(v,  lU  n 

H  (S51/2) 


+  II  (v 


f-1  r0H°(S„  ) 

o(/2 


From  the  induction  hypothesis, 

,2 


|V(-1V'2(s  ,  ’  "v 

%  <S5,/2) 


**'-,<v.wi„«s4i/2, 


*'ir 


<  00. 


Since  r  <  1  on  S  „  we  see  that  ||vj  , 

5l/2  1  n\s^  ) 


C  ||  v 


£-1' 


5 ;  /2 


i«2. 2 

p  o ,  /2  ‘ 


and 


hence  that  II  v 


,1, 


»ii  .  <  oo.  Now,  since  v  e  H  (S  )  and  because  of  (4.3) 

H1  ( S„  )  1  5‘ 

O  [  /  tL 


0 


we  can  apply  Lemma  2.3  to  get,  for  0  <  f3  <  1  satisfying  (3.  <  13, 


<  c. 


1-2  2 
A  ||  r  i  r  u ) 


l'1' £  (S  )  +  ^  V£^  l  1 

r*VV 


,5^)|a||i 


<  C  max(  1 , 5'.J) A||r^  2(r2u)  J 

1  11  r  H  (S, 


+  v 


i  1 

K  (S5i'S51/2) 
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i.iA.n.t, 


=  Cjjxilr*  2(r2u) 


e!l  0  -  II  v  i 

r  H  (S5i)  ^(S5i-S5i/2)- 


Then  we  note  that 


(4.6) 


r£_2(r2u)  =  V2£v£-l  +  *U~1)v£-2’ 

r 


(4.7) 


.  1-2,  2  , 

I r  (r  u) 


r  o 

r  HU(S-  ) 

°i 


<  ||r£  2(r2u)  J  +  II r£  2(r2u)  J  Q 

r  HU(Sa  -S.  ._)  r  H(S  ) 

O  i  O  [/  Oj/  c. 


(4.8) 


iv  H  =  ||(l-£)v  +r(v  )  || 

H  (S5l/2)  H  (S61/2) 


5  C«-DI!v  H  +lr(v  )  || 

h°(s5i/2)  h°(s6i/2; 


S  £||  v 


£- 1  1 

*  <Sa,/2  > 


S  £||  v 


£-l"^,2 

*0  S51/2) 


Now,  using  (4.3),  (4.6),  (4.7),  (4.8),  and  the  induction  hypothesis,  we  obtain 


1-2  2 
|r  “(r  u) 


£“  o 

r  HU(S, 


<  || 2(r2u)  .||  +  l|r£  2(r2u)  J  Q 

r  H  (sarsa,/2)  r  H  (s5l/2: 


S  II  v  |  +  2£||  v  || 

H  (S6i-S6i/2)  H  (S61-S61/2) 


♦  £(£-l)||v 


£-2  0 

H  ls5,-sa,/2> 


>  V  V  V  V  V  V  V  V  V  V  V  V  V  V 

a.  -N  »V  V 


v: 


Note  that 
(4.3),  (4. 

(4. 10) 


+  «V  o  +2£||v  || 

h°(s5i/2)  h°(s5i/2) 


+  £(£-l)||v 


P-2  n 

H  (sa,/2> 


<3£||v  II  +  £(£-1 )  ||  v.  II 

^  (S5l/2}  ^  (S61/2) 


+  II  v  ||  +  2£||v  || 

*  (S8i-S51/2)  '  K  ( W2} 


+£(£-l)||v 


P-?"  1 

*  <W 2> 


£-1  £-2  £  £-1  £-2 
<£![C(3d  i+a  )  +  CQ(dQ+2dQ  +d*  )]. 


/ 

/ 

i 

3“ 


1-2  t-2 

if  £  =  1,  the  terms  d  and  dQ  are  not  present.  Combining 
5),  and  (4.9)  we  get 


£"^2,2 

%  ( S«,/2' 


£-1  t-2 

<  £!{CC1A(3d  +d  ) 


*cociMdo*2do"*do'2)*cocidol 


=  £!d£|c 


,3C  A  CjA-» 


+  CqC^A+1) 


+  2CoClA 


+  cocix 


'V"2 


d2} 


<  £!  d  <C 


,3C  A  C  A-» 


C_C. ( A+l )d„  2C„C. A  C„C„A> 


•T  V*  -rv-:r> 


<  £!  d 


Tc  c' 
[2  +  2_ 


=  C£!d  ; 


thus  (4.4),  for  k  =  i,  holds  with  C  and  d  as  selected  above.  This 
completes  the  proof  of (4. 4)  for  all  k  >  0. 

Next  we  show  that  there  are  constants  C  and  d  such  that 


(4.  11) 


,  cti-2~a  ,  --k, 

ilr  X  mi  .  <  Cd  k!  , 

5i/2; 


for  all  a  with 


I  a  I  =  k+2  >  2  and  <  2, 


provided  (3^  <  3-  We,  in  fact,  show  that  (4.11)  holds  with  C  =  3C  and 
max(6,d),  where  C  and  d  are  the  constants  in  (4.4).  Let  k  >  0.  Then 
there  are  three  choices  for  a  with  |a|  =  k  +  2  and  <  2,  namely  a  = 

( a  1 ,  <*2 )  with  =  k  and  -  2,  =  k+1  and  <*2  =  1,  and  0^  =  k  +  2 

and  =  0.  Thus  (4.11)  is  equivalent  to  the  following  three  inequalities 


(4.  12a) 


(4.  12b) 


(4. 12c) 


Irk"2u  k  zl|  }<Cdkk!, 

r  9  V  5i/2J 


<s,  „)  s  E5kk!’ 

p  0[/Z 


!rku  k*2|l  <Cdkk!,  k  >  0. 

T6|/2 


Recall  that  (4.4)  states  that 


(4.13)  ||  v 


11 2  m  k  ..2  n  (3.  k  1  1, 2 

k"„2-2(s  “  |r  urk  1  *  lr  I-  V'rr''  0.„  . 

^  S6,/2  S5 ( /2  H  S5,/2' 

3-1,  k  ,  m2 

*!lr  <r  V’rel.O.-  . 


ai 


(4.12a)  follows  immediately  from  (4.13).  Next,  consider  (4.12b).  For  k  = 
it  follows  directly  from  (4.13).  Suppose  now  that  (4.12b)  holds  for 
0,1,..., k-1.  Then,  using  (4.13),  we  have 


,  k-1  ,  ,,  13-1.  k  .  ,  £+k-2 

lr  ^'»"VV'lr  urk  re  Ve«„0(V2, 

\c  — lc-1 

<  Cd  k!  +  Cd  k! 
k  -k-1 

=  Ck! (d  +3d  ) 


<  2Ck!dk 


<  Cdkk! , 


which  shows,  by  induction,  that  (4.12b)  holds  for  all  k.  Finally  consider 
(4.12c).  For  k  =  0,  this  follows  directly  from  (4.13).  Suppose  now  that 
(4.12c)  holds  for  0,1,..., k-1.  Then,  again  using  (4.13),  we  have 


lrku  k*2||  ,  }  =  llrArSi  k) 

/T  Si/2 


-  2kr0+k_1u  k+1  -  k(k-l  )rP+k_2u  k|| 

H  (S61/2) 

V-  --k-1  — k-P 

<  Cd  k!  +  2Cd  k!  + Cd  k! 


,  ,  ,k  -k-1  „-:k-2, 
=  Ck! (d  +6d  +3d  ) 


<  (3Cik!dk 


=  Cdkk! , 


which  shows  that  (4.12c)  holds  for  all  k.  This  completes  the  proof  of 


Using  the  identity 


A( raiu 


r«iea2-2 


)  =  rai  2(r2Au  ) 

ea 2-2  rai 


.  ai-2,  2  . 

=  "Ar  (r  u  ) 

0  2-2  rai 


which  is  valid  for  ^  2,  we  get 


(4. 14) 


ai-2 

r  u 


raiea2 


-(rttlu 


a,  a2-2  2 

r  6  *  r 


“1/  ai  \ 

-  r  ( r  ‘u  „) 

ai„a2-2  r 
r  Le  2 


cti~2 


-  Ar  (r  u  ) 

q“2-2 


Note  that  the  number  of  0-differentiations  on  the  right  side  of  (4.14)  is  two 
less  than  on  the  left.  Thus  repeated  application  of  (4.14)  reduces  the  number 
of  0-differentiations  to  <  2,  so  that  (4.11)  can  be  applied.  In  this  way 
one  shows  that  there  exists  C  and  d  so  that 


(4.  15) 


|rai  22)au|| 


==k 


<  Cd  k! ,  for  |a|  =  k+2  >  2, 


0(S 


5t/2; 


provided  0^  <  0. 


Now  u  €  $2(S.  )  if  and  only  if  u  e  H2,2(S  )  and 

p  o  \/  d  p  0\/d 


(4. 16) 


\rl3+k  2Dau||  <  Cdk  2(k-2)! 

"  lS5,/2 


I  a |  =  k,  k  =  2.3. . . 


for  some  C  and  d  independent  of  k.  (4.4),  with  k  =  0,  together  with 


2  2 
rc'  »  c 


Lemma  2. 1,  shows  that  u  e  H''  (ST  j.  It  thus  remains  to  prove  (4. 16).  L 

p  O  \  /  d 


k  >  2  be  fixed.  By  Lemma  2.2, 


i  /3+k-2  a 
r  D  u 


<  Cd  (k-2) ! .for  |a|  =  k 


H  (S61/2) 


if  and  only  if 


(4.  17) 


k' -2-a'  , 

I r  2Da  ul|  <  Cd  (k'-2)!  , 

V  a  1/2 


for  2  <  |a' |  =  k'  <  k. 


But  from  (4.15)  we  have 


k'-2-a'  , 

I  r  V  ully  ,  *  Cd*  ~^(k'-2)! 

8l  6,/ 2J 


for  2  <  |a' |  =  k'  <  k, 

which  implies  (4.17).  Thus  (4.16)  holds  for  some  C  and  d,  and  we  conclude 

2  0 
that  u  e  2_(S_  ._)  for  any  0  <  ft  <  1  satisfying  1  ~  n/u  =  (3  <  13. 

13  0\  /  cL  1  i 


In  this  way  we  prove  that  u  e  (S„  .„)  for  fi  <  ft.  in 

p  i  O  i  /  d  1  1 


the  ne  i  p.hLcr- 


hood  of  each  vertex  A. .  Combining  this  result  with  Lemma  4.1  we  arrive  at 
the  proof  of  Theorem  3.  1  in  the  case  that  L  =  ~A  and  tl  is  a  straight 
polygon.  If  L  is  the  general  operator  in  (3.2)  and  Q  is  a  straight 
polygon,  the  result  is  obtained  using  the  techniques  in  Section  3.2  of  [5], 
Finally,  for  the  general  case  in  which  n  is  a  Lipschitz  curved  polygon  wo 
employ  the  arguments  used  in  the  proof  of  Theorem  3.4  of  [5]  to  get  the 

2 

desired  result.  More  specifically,  by  that  technique  we  show  that  u  6 
Then,  using  Lemma  2. 4, we  get  Theorem  3.1  in  its  full  generality. 

Remark  4.1.  A  careful  analysis  of  the  proof  of  Theorem  3.1  allows  one  to 
assess  the  dependence  on  A  of  the  values  at  the  constants  C  and  d  (cf. 
(4.16)),  which  are  related  to  the  smoothness  of  the  eigenfunctions.  One  car. 
see,  for  example,  that  d  =  dA + k,  where  d  >  1  and  k  are  independent  of 


A.  This  suggests  that  the  higher  eigenfunctions  are  less  smooth  than  the 


lower  ones.  Although  this  is,  in  general,  correct,  we  shall  see  in  Section  6 
that  there  are  important  exceptions.  In  the  example  treated  there  we  will  see 
that  of  the  first  three  eigenfunctions,  the  first  is  the  roughest  while  the 
third  is  the  smoothest. 


5.  Basic  Results  on  Eigenvalue  and  Eigenfunction  Approximation 


As  we  have  seen  in  Section  3,  the  eigenvalue  problems  we  are  considering 


have  the  form 


(5.  1) 


Seek  A  e  R,  0  *  u  e  Hq(Q)  satisfying 
a(u,v)  =  Ab(u,v)  V  v  e  H^CQ), 


where  the  bilinear  forms  a  and  b  are  defined  in  (3.7)  and  (3.8). 

We  are  interested  in  approximating  the  eigenpairs  (A,u)  of  (5.1),  and 
toward  this  end  we  select  a  family  of  finite  dimensional  subspaces  S(p,k)  of 
Hq(Q),  indexed  by  the  parameters  p  =  1,2,...  and  k  =  0,1,...,  satisfying 


(5.2) 


inf  ||u-*|| 
sSlp.k)  H 


>0  as  p — >»,  uniformly  in  k, 


for  each  u  e  HQ(n) 


(the  reason  for  employing  two  parameters  will  be  made  clear  later).  Then  we 
consider  the  finite  dimensional  eigenvalue  problem 


(5.3) 


Seek  A(p,k)  €  R,  0  *  u(p,k)  e  §(p,k)  satisfying 
a(u(p,k),v)  =  A(p, k)b(u(p, k) , v)  V  v  e  §(p,k). 


(5.3)  has  a  sequence  of  eigenvalues 


0  <  A  (p, k)  <  ...  i  AN(p, k) 


and  corresponding  eigenfunctions 


sat isfying 


u1  (p,  k) . uN(p,  k) 


a(u.  (p,  k) ,  u  .(p,  k) )  =  5.  ., 
i  J  ij 


where  N  =  N ( p , k )  =  dim  §(p,k).  The  eigenpairs  ( A  (p, k) , u  ( p, k ) )  of  (5.2) 
are  viewed  as  approximations  to  the  eigenpairs  (A.,u.)  of  (5.1).  It  is 


well-known  (see,  e.g.,  [7])  that 

(5.4)  A.  <  A.(p,k),  V  p, k,  i  =  1 . N, 

and,  as  a  consequence  of  (5.2)  (see,  e.g.,  [7]),  that 

lim  A.(p,k)  =  A.,  uniformly  in  k,  for  each  i. 

p-*oo 

We  now  present  the  error  estimates  that  describe  the  quality  of  the 
approximation  (A^(p, k) , u  (p, k) ) .  For  the  sake  of  simplicity,  we  state  these 
results  in  the  case  in  which  the  eigenvalues  A^  of  (5.1)  are  simple.  For  a 
complete  treatment  of  the  general  case  we  refer  the  reader  to  [7,8,9,10,11]. 
Let 

l|u|iE  =  Va(u,  u) 

and 

||u!!b  =  \/b(  u,  u) . 

Then  we  let 

(5.5)  c^p.k)  =  inf  Hu^zllg. 

*€§(p, k) 

(5.6)  i>(p,k)  =  sup  inf  li  Tf-^li  E» 

feH^(n)  xe§(p,k) 

!|fHE=l 

and 


5.7) 


T7(p,k)  =  sup 


inf 


ITf-xl1 


f€L2(fl)  *e§(p,k) 


E’ 


!f"0=1 


where  T  :  (-^(fl) — »H  (fi)  is  the  solution  operator  corresponding  to  the 


differential  operator  L  introduced  in  (3.1),  i.e.,  T  is  defined  by 


The  error  estimates  are  given  in  terms  of  the  quantities  e^(p,k),  i/(p,k), 
and  T)(p,k).  We  note  that  all  of  these  quantities  approach  zero  as  p — >oo, 
uniformly  in  k. 

Theorem  5.1.  Suppose  A^  is  simple.  Then  there  is  a  constant  such  that 

(A  (p, k)-A. )/A 

(5.9)  1  -  d.T)  (p,  k)  <  — - — 5 - - - -  <  1  +d.w(p,k), 

1  cf (p, k)  1 

for  all  p  and  k, 

and 

II u  (p,  k)-u.  ||  /|| u  || 

(5.10)  1  s  ~  £  1  *V(P.k>. 

for  all  p  and  k. 

This  result  shows  that  for  k  fixed  and  p  large,  the  eigenvalue  error 

(A. (p,k)  -  A. )/A. 

is  nearly  equal  to 

e^(p, k) , 

and  the  eigenfunction  error 

||u.  ( p,  k)  -  u^lg/Hu^lg 

is  nearly  equal  to 

c. (p, k) . 
l 

One  of  the  goals  of  this  paper  is  to  make  a  computational  study  of  these  error 
assessments  in  the  practical  range  of  the  parameters  p  and  k. 


By  the  usual  duality  argument  we  have 


(5. 11) 


f|ui  -  E(p,  k)u.||b  <  Cr)(p,  k)||ui  -  E(p,k)uiHE> 

where  E(p,k)u  is  the  a-orthogonal  projection  of  onto  §(p,k),  showing 

that  ||u.  -  E(p,k)u  ||  =  o(  ||u^  -  E(p,  k)u.  ||  ^.) .  A  parallel  result  holds  for  the 

eigenfunction  error  u.  -u.(p,k).  We  state  this  in 

Theorem  5.2.  Suppose  A.  is  simple.  Then  there  is  a  constant  C.  such  that 


(5. 12) 


|u.  —  u.  ( p,  k )  ||  b  <  C.T)(p,k)||u.  -  u.Cp.kJUg.. 


This  result  is  proved  in  [11]. 


6.  Numerical  Results  and  Their  Analysis 

In  this  section  we  consider  a  model  problem,  present  numerical  results 
for  it,  and  analyze  these  numerical  results  in  light  of  the  theoretical 
results  outlined  in  Section  5.  We  will  consider  the  behavior  of  both  the  p 
and  h-p  versions  of  the  finite  element  method.  The  specific  aim  is  to 
assess  the  reliability  of  the  theoretical  results  -  which  are  of  an  asymptot 
nature  -  as  a  guide  to  practical  computations  -  which  may  take  place  in  the 
pre-asymptot ic  phase  -  and  to  look  for  characteristic  features  of  the  numeri 
cal  results  which  are  not  completely  explained  by  known  theoretical  results. 

The  Model  Problem 

Let  £2  be  the  L-shaped  domain  shown  in  Figure  6. 1  and  consider  the 
eigenvalue  problem 

f-Au  =  Au,  x  €  £2 

(6‘  U  |u  =  0,  x  e  T  =  3 £2. 


Figure  6.1.  The  Domain  £2 


n  is  a  Lipschitz  straight  polygon,  as  described  in  Section  2,  and  (6.1)  is  an 
eigenvalue  problem  of  the  type  described  in  Section  3.  We  are  taking  =  r. 
The  bilinear  forms  a  and  b  are  given  by 


a(u, v)  = 


fau 

a  v  au 

3v  ' 

K 

n 

SX2 J 

(6.3) 


b(u,v)  =  uv  dx. 


and  the  norms  ||u||^.  and  ||u||^  are  given  by 


(6.4) 


(6.5) 


!!  u||  b  =  u2  dx 


Thus  || u || is  the  usual  energy  norm  in  ( C2 )  and  ||u|lb  the  usual  L^-norm. 


0<A1<A2<...  /co 


U 1 ’ U2 ’ 


be  the  eigenvalues  and  eigenvectors  of  (6.1)  and  suppose 


a‘Vuj)  =  6i.r 


The  model  problem  (6.1)  was  selected  as  a  typical  example  of  a  problem 


a  domain  with  nonsmooth,  but  piecewise  analytic,  boundary.  It  will  also  be  of 


interest  in  that  it  has  some  symmetries.  (6.1)  is  a  reasonable  representation 
of  a  large  class  of  engineering  problems. 

The  Regularity  of  the  Eigenfunctions 

On  the  basis  of  Theorem  3. 1  we  can  make  an  assessment  of  the  regularity 

of  the  eigenfunctions  of  (6.1).  It  is  easy  to  see  that  the  eigenfunctions  of 

(6.1)  are  analytic  at  the  vertices  A.,  where  w.  =  n/2,  j  =  2 . 6,  but 

J  J 

3 

that  at  A^,  where  =  ~n,  the  eigenfunctions  are  singular.  From  Theorem 
3.1,  more  specifically  from  the  results  of  Section  4,  we  see  that 


(6.6) 

u.  = 

l 

where 

0  = 

(3r. 

••’V 

)  with 

0  <  %■■■ 

.  .  ,/3_  <  1  arbitrary  (in  fact,  we  can 
b 

take 

02... 

••*6 

=  0) 

and  1 

-  —  =  i  < 

3 

We  can  also  use  the  more  standard  regularity  theory  for  problems  in 
domains  with  corners  (see,  e.g. ,  [12])  to  analyze  the  singularities  of  the 
eigenfunctions.  We  have,  for  example, 


p/o 

(6.7)  u  =  (CjP  sin  2/3  e)^  +  i>r 

00  2 
where  i//  is  C  cut-off  function  based  at  the  point  A^,  €  H  (fi),  (r,0) 

are  the  polar  coordinates  depicted  in  Figure  6.1,  and  *  0.  Obviously 

a  5/3  5/3 

u,  e  H  (fi),  with  a  <  1+2/3,  but  u,  «  H  (Q).  Nevertheless  u.  e  2L  (Q). 
1  1  1  2, 00 

(Here  Ha(n)  is  the  fractional  Sobolev  space  and  B*  (Q)  is  the  Besov 

C*  ,  00 

space;  see  [13]).  Not  all  the  eigenfunctions  have  the  form  (6.7).  In  fact, 
symmetry  considerations  show  that  is  antisymmetric  with  respect  to  the 

line  x  =  -y  and  hence 

4/3 

(6.8)  u2  ~  ^C2r  Sin  + 

where  v ^  €  H^(Q)  and  *  0,  so  that  u^  e  Ha(fi),  for  any  a  <  1  +4/3, 


28 


and  u„  e  B^+4/3(Q).  Furthermore, 

2  2,  oo 

axes  and  thus 


u^  is  antisymmetric  with  respect  to  both 


(6.9) 


,2.1/2  1  , 

u_  =  (-)  -  sin  rrx,  sin  mx  • 

J  J  71  1  2 


u^  is  therefore  analytic  on  fi. 

We  see  here  a  typical  phenomena  for  eigenvalue  problems  arising  in 
engineering,  namely,  because  of  the  presence  of  various  symmetries,  higher 
eigenf unct ions  may  be  smoother  and  hence  lead  to  more  accurate  numerical 
approximations  than  the  lower  ones.  See  also  Remark  4.1. 


The  Finite  Element  Spaces 


We  consider  the  meshes  shown  in  Figure  6.2.  These  are  typical  meshes  for 
the  p  and  h-p  versions  of  the  finite  element  method  for  the  approximation 
of  (6.1)  since  the  leading  singulari t ies  of  the  eigenfunctions  are  located  at 
vertex  (for  more,  see  [2]  and  [14]).  The  meshes  are  characterized  by  the 

number  k  of  layers. 
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v.;.; 


v  j*  v*  •  * 


v.v.v. 


By  S(p,k)  we  denote  the  finite  element  space  consisting  in  polynomials 


of  degree  p  on  a  mesh  with  k  layers.  The  functions  on  the  trapezoidal 
elements  are  defined,  as  usual,  as  pull-back  polynomials  of  degree  p  based 
on  the  standard  bilinear  map.  The  elements  are  of  serendipity  type  Q'  (see 
(15]).  The  number  of  degrees  of  freedom  of  one  (unconstrained)  quadrilateral 
element  is 


4p,  for  p  <  4 
4p  *  lp-^(p~3). 


for  p  >  4. 


By  §(p,  k)  we  denote  the  finite  element  space  constrained  by  zero  on  F, 
i.e.,  §(p,k)  =  S(p,k)  n  H^(n).  The  number  of  degrees  of  freedom  is  defined 
to  be  the  dimension  of  §(p,k)  and  will  be  denoted  by  N(p,k).  For  large  p 
and  k  we  obviously  have 


N(p,  k)  *  (k+1  )p2. 

For  small  p  and  k  this  asymptotic  formula  is,  of  course,  not  very  accu¬ 
rate.  Table  6.1  gives  N(p,k)  for  various  combinations  of  p  and  k  under 

2 

consideration.  The  ratio  N(p, k)/(k+l)p  ranges  from  .44  to  3.75  over 
our  range  of  p  and  k. 


1 

2 

3 

4 

5 

6 

8 

0 

0 

2 

4 

9 

17 

28 

42 

59 

1 

5 

18 

31 

53 

84 

124 

173 

231 

2 

10 

34 

58 

97 

151 

220 

304 

403 

3 

15 

50 

85 

141 

218 

316 

435 

_ 

575 

Table  6.1.  Number  of  Degrees  of  Freedom 


We  recall  that  in  the  p  version,  accuracy  is  achieved  by  letting  p — >c 
while  k  is  held  fixed  (i.e.,  the  mesh  is  held  fixed).  With  the  h-p 
version,  we  simultaneously  increase  p  and  k.  In  the  classical  h  version 
accuracy  is  achieved  by  refining  the  mesh  while  keeping  p  fixed.  We  will 
concentrate  on  the  p  and  h-p  versions.  We  will  not  study  the  h  version 
nor  assess  the  question  of  how  the  performance  of  elements  of  type  Q' 
(serendipity  type),  which  we  are  using,  compares  with  that  of  the  full  tensor 
product  elements  (of  type  Q)  or  with  that  of  triangular  elements.  Neverthe¬ 
less,  we  will  consider  the  effect  of  distortion  of  elements  for  the  p  and 
h-p  versions.  For  k  =  0,  none  of  the  elements  are  distorted,  while  for 
k  >  0  many  elements  are  trapezoidal,  and  hence,  distorted. 


In  (5.5)  we  introduced  the  quantity 

c.(p.k)  =  inf  ;|u.-*ilE; 

*€§(p, k) 

clearly 

c.  (p,  k)  =  || u .  —  E(p,  k)u.  || _ 

1  1  1  C 

||u.-E(p,k)u.  || £ 


C 

where  E(p,k)u  is  the  a-orthogonal  projection  of  u^  onto  S(p,k)  (recai 
that  ||u^||E  =  D-  We  have  seen  in  Section  5  that  the  accuracy  of  the  finite 
element  approximation  of  the  eigenpairs  is  determined  mainly  by  e  (p,k).  We 
now  present  some  theoretical  and  numerical  results  on  the  size  of  e  (p,k).  W 

first  note  that  (5.2)  holds  for  our  choice  of  §(p,k).  To  see  this  we  note 

00  1 
that  the  C  functions  with  compact  support  in  n  are  dense  in  H  .  Hence 

condition  (5.2)  is  equivalent  to  the  approximation  of  smooth  functions  by  the 


p  version  and  the  results  of  [16]  lead  to  the  desired  conclusion. 


From  the  approximation  results  in  [16]  and  the  regularity  results  (6.7, 


(6.9)  we  obtain 


(6. 10a) 


(6. 10b) 


(6. 10c) 


Cjtp, k)  <  CjCklp 


e2(p, k)  <  C2(k)p 


c  (p.k)  <  ( k ) p  ,  for  any  c  >  0. 


These  are  estimates  for  the  p-version  of  the  finite  element  method  on  a  fixe 
mesh;  (6.10a)  and  (6.10b)  are  optimal,  but  (6. 10c)  is  not.  In  fact,  c  (p,k) 
decreases  exponentially  in  p.  The  estimates  (6.10)  (including  the  ref ir.emer 
of  (6.10c)  just  mentioned)  can  obviously  be  written  in  terms  of  N(p,k).  If 
this  is  done  e^(p,k)  and  e2(p,k)  are  seen  to  decrease  algebraically  and 
c  (p.k)  exponentially  in  N(p,k). 

If  p  2.  1  and  k  >  0  are  related  by 

( k+ 1 )  <  p  <  a2 ( k+ 1 ) , 


where  0  <  a ^  <  a^,  then  (see  [5]  and  [14] 


(6.11)  c.  (p, k )  <  Co  I(NfP'k11 


,  for  some  C  >  0  and  y  '•  0. 


j  and  C,  in  general,  depend  on  i,a  and  the  domain,  but  are  indupor 

2 

dent  of  N.  (6.11),  which  follows  from  the  fact  that  u.  e  2^(0)  for  every 

i  (cf.  (6.6)),  is  the  basic  estimate  for  the  h-p  version. 

Figure  6.3  depicts  the  relative  error  c  (p.k)  for  various  k  as  a 

-4/3 

function  of  p  in  double  logarithmic  scale.  The  rate  p  is  shown  as  a 

slope  in  the  figure  (cf.  (6. 10a)).  We  see  the  typical  behavior  of  the  p 


version  of  the  finite  element  method. 


."V  V  V  K.TF.”*.’ V". V  V 


The  reverse  S-curve  shown  in  Figure  6.4  characterizes  the  performance  of 


the  p  version.  It  has  two  parts,  the  pre-asymptot ic  phase  where  the  curve 


turns  down  and  essentially  gives  the  behavior  of  the  h-p  version  and  the 


asymptotic  phase  where  the  curve  is  a  straight  line  (with  slope  -4/3  in  cur 


case).  As  k  increases  the  asymptotic  phase  shifts  toward  higher  p  (see 


Figure  6.3). 


To  understand  this  behavior  we  note  that  the  asymptotic  phase  begins  with 


that  p  at  which  the  error  due  to  the  elements  with  vertices  at  (where 


the  leading  singularity  occurs)  starts  to  dominate  the  total  error.  Because 


u.  €  STTn),  the  error  decreases  exponentially  in  the  pre-asymptot ic  phase 
i  P 


during  which  it  is  not  dominated  by  the  elements  at  A^.  We  further  note  that 


for  k  low  (i.e.,  for  an  unrefined  mesh)  the  singularity  greatly  affects  the 


error  and  thus  its  influence  shows  up  already  for  low  p,  whereas  for  k 


high  (i.e.,  for  a  refined  mesh)  the  singularity  has  been  resolved  by  the  mesh, 


hence  has  less  effect  on  the  total  error,  and  its  influence  shows  up  only  for 


high  p.  Roughly  speaking,  the  mesh  is  not  properly  refined,  for  the  desired 


accuracy,  if  we  are  in  the  asymptotic  phase  of  the  p  version  (and  the  rate 


of  convergence  is  algebraic).  For  more  details  see  [2] 


Figure  6.5  presents  the  relative  error  e^Cp.k)  for  various  k  as  a 


function  of  p.  Here  the  slope  in  the  asymptotic  phase  of  the  p  version  is 


-8/3  (cf.  (6.10b)).  For  k  =  0  we  are  in  the  asymptotic  range  for  p  >  6, 


but  for  k  =  2  and  3  the  asymptotic  phase  is  not  reached  for  the  values  of 


p  we  are  considering. 


V  V.V.  r'*  ’  \.r' 1 


Figure  6.3  shows  an  interesting  phenomena  occurring  in  practice,  namely, 
that  in  certain  circumstances  refinement  can  lead  to  an  increase  in  the  error. 
We  note  that  for  fixed  k,  the  inclusion  §(k,p+l)  o  §(k,p)  is  valid  and 
this  leads  to  a  decrease  in  error  with  increasing  p,  while  for  p  fixed  we 
do  not  in  general  have  §(k+l,p)  o  §(k,p)  so  that  the  error  is  not  guaranteed 
to  decrease  with  increasing  k.  For  low  p,  the  major  error  occurs  in  the 
largest  elements  and  in  these  elements  distortion  deteriorates  approximate  1 i ty 
in  H 1 ( f2 ) .  A  priori,  it  is  virtually  impossible  to  predict  this  behavior.  In 
general,  a  conservative  strategy  for  mesh  refinement  is  probably  advisable. 

If  one  is  interested  only  in  very  low  accuracies,  in  certain  situations  an 
unrefined  mesh  with  undistorted  elements  could  give  better  results  than  a 
refined  mesh  with  distorted  elements.  Comparing  Figures  6.3  and  6.5  we  see 
the  interesting  feature  that  for  p  =  2  the  refined  mesh  gives  a  better 
result  for  the  second  eigenfunction,  while  for  the  first  eigenfunctions  the 
unrefined  mesh  gives  the  better  result.  Because  we  usually  compute  several 
eigenvalues  with  one  mesh,  we  see  that  the  conservative  refinement  strategy 
mentioned  above  is  advisable.  (We  remark  that  for  triangular  meshes  the 
situation  is  different  since  in  this  case  the  spaces  are  nested. ) 

Figure  6.6  shows  the  behavior  of  e^Cp.k).  Since  u^  is  analytic,  the 
exponential  rate  occurs  for  all  p  (cf.  (6.10c)  and  the  following  discus¬ 
sion)  and  essentially  no  difference  can  be  seen  in  e^tp.k)  as  a  function  of 
k  >  1.  It  is  thus  best  to  use  an  unrefined  mesh. 


Since  for  k  =  0,  is  symmetric  in  every  element  of  the  mesh,  the  error  for 

odd  p+  1  is  the  same  as  that  for  even  p.  We  note,  however,  that  in 
general  such  symmetries  are  difficult  to  predict  a  priori. 

Remark  6.1.  The  numerical  results  we  present  require  the  values  of  the  exact 
eigenfunctions,  which  are  not  explicitly  known.  We  obtained  values  for  these 
eigenfunctions  by  a  careful  extrapolation  procedure  which  we  believe  leads  to 
reliable  results. 

Accuracy  of  the  Approximate  Eigenvalues 


Consider  the  finite  dimensional  eigenvalue  problem  (5.3),  where  the 
spaces  §(p, k)  are  as  described  in  the  Subsection  The  Finite  Element  Spaces. 
With  this  choice  for  §(p,k),  (5.3)  defines  the  finite  element  method  for  the 

approximation  of  the  eigenpairs  of  (6.1).  As  in  Section  5,  denote  by 
(Ai (p,  k) , u. (p,  k) )  the  eigenpairs  of  (5.3).  These  (A. (p, k) , u. (p, k) )  are 
called  the  (finite  element)  approximate  eigenpairs.  As  noted  in  Section  5, 

A.  <  A. (p, k) ,  V  p, k,  i  =  l,...,N(p,k)  =  dim  §(p,k), 


lim  A.(p,k)  =  A.,  uniformly  in  k,  for  each  i. 

p-»oo 

We  are  interested  here  in  the  accuracy  of  the  approximation  A^(p,k)  - 
A.;  specifically  we  are  interested  in  the  comparison  of  (A^ (p, k)-A^ )/A^  and 
e?(p,k)  (cf.  (5.9)).  Thus  we  define 


(6. 12) 


C  (p,k)  = 


(Ai(p,k)-A.)/A. 

e2(p, k) 
i 


From  Theorem  5.  1  and  the  facts  that  in  our  case  r(p,k)  <  0(p  /  )  and 
Tj(p,k)  <  0(p  ^),  uniformly  in  k,  (cf.  [16]),  we  have 


-2  -4/3 

l-0(p  )  <  C,(p,k)  <  l+0(p  ),  uniformly  in  k. 
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Table  6.2  gives  the  values  of  Ch(p,k)  for  various  values  of  p  and  k 


gjl 

IS 

0 

1 

2 

3 

1 

— 

1. 322 

1.349 

1. 346 

0.912 

0.  988 

0.  978 

0.  978 

B 

0.967 

0.947 

0.932 

0.  932 

B 

0.962 

0.945 

0.929 

0.927 

1 

5 

0.977 

0.  991 

0.  965 

0.  954 

6 

0.984 

0.  995 

0.973 

0.  362 

7 

0.992 

0.998 

0.  997 

0.  992 

8 

0.995 

0.999 

0.  998 

0.  996 

1 

- 

1.795 

1.777 

1.777 

2 

0.867 

0.  919 

0.  921 

0.978 

3 

0.787 

0.892 

0.893 

0.932 

4 

0.780 

0.  886 

0.886 

0.927 

d 

5 

0.871 

0.932 

0.931 

0.934 

6 

0.948 

0.932 

0.929 

0.962 

7 

0.966 

0.963 

0.962 

0.991 

8 

0.975 

0.974 

0.975 

0.995 

1 

- 

2.692 

2.740 

2.740 

2 

- 

1.  126 

1.  128 

1.  128 

3 

1.  150 

1.  121 

1.  124 

1.  124 

4 

0.831 

0.902 

0.  902 

0.901 

J 

5 

0.831 

0.899 

0.  899 

0.  898 

6 

0.905 

0.904 

0.905 

0.  904 

7 

0.905 

0.928 

0.928 

0.928 

8 

0.942 

0.  946 

0.946 

0.946 

Table  6.2.  The  Values  of  C^(p,k). 

We  see  that  if  e^(p,k)  <  5%  (cf.  Figures  6. 3, 6. 5, 6. 6) ,  then  C.(p.k) 
and  Cj(p,k)  <  1,  and  furthermore,  C^(p,k)— >1  as  p — >oo. 


The  fact  that  C.(p,k)  <  1  for  most  values  of  p  and  k  can  be 


explained  in  part  as  follows.  Let 


E(p,k)u  =  y  a(u, u  .(p, k) )u  (p, k) , 
l  t  J  J 

j=l 

be  the  a-orthogonal  projection  onto  §(p,k)  and  let 

N 

E^lp.klu  =  ^a(u,  Uj(p,  k)  )Uj(p,k). 


J=i 


Then  we  have 


(6. 13a) 


(A. (p.k)-A. )/A. 

c^(p, k) 


f||ui-E.(p,k)uil||  ||ui-Ei(p,k)u.||^'j 


1 


[a i  || E j,  ( p ,  k ) u i  || ^  ||Ei  (p.  k)uit,b 


2  (  2. 

c^P. 


1  -  A 


lUi-EjIp.klUi^' 


i||ui-Ei(p,k)u.||gJ 


lUj-E.  (p,  k)Uj  ||  j. 


1 


Ail|ui-E(p.k)uill^  ||E.  (p.  k)ut  ||b 


Note  that  for  i  =  1,  (6.13a)  simplifies  to 


(6. 13b) 


(AJ (p, ki-Aj ) /A  ^ 

e^( p, k ) 


1  -  A, 


|u1-E(p,  k)u1!lb 


j|u1“E(p,  kju^, E 


||-x{A1l|E(p>k)u1l|^}"1 


Table  6.3  shows  the  values  of  the  right  side  of  (6.13b)  for  k  =  0 


Table  6.3.  The  Values  of  the  Right  Side  of  (6.13b)  (for  A ^ ) . 

We  note  two  features: 

•  The  values  in  Table  6.3  are  very  accurate  upper  bounds  for  the  correspond¬ 
ing  upper  bounds  in  Table  6.2.  This  indicates  that  the  inequality  in 
(6.13b)  is  a  near  equality. 

•  There  are  two  factors  on  the  right  side  of  (6.13b);  the  first  is  less  than 


1  and  the  second  may  be  greater  than  1.  Since  the  values  in  Table  6.3 


are  mostly  less  than  1  we  see  that  the  first  factor  on  the  right  side  of 


(6. 13b)  is  having  a  greater  influence  than  the  second  factor.  Thus  the 


cancellation  in  the  first  factor  is  causing  the  right  side  of  (6.13b)  to  be 


less  than  1  and 


|u^-E(p, klUjU^ 


1l|u1-E(p,k)u1ll| 


plays  the  major  role  in  determining  the  degree  to  which  it  is  less  than  1. 


Although  we  don’t  have  computed  values  for  the  right  side  of  (6.13a), 


Table  6.4  shows  the  values  of 


||u  -E(p,k)u  || b]  2 

l-A.—i - ||E(p.k)u.O 

||ui-E(p,k)u.  || £J 


i 

i 
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for  i  =  3.  We  believe  this  expression  is  very  close  to  the  right  side  of 


(6.13a)  since  E^(p,k)u^  -  E(p,k)u^.  We  again  see  the  two  features  noted 


above. 


\\ 

p  \ 

0 

1 

1 

- 

- 

2 

- 

6.52 

3 

1. 28 

4 

0.831 

0.  932 

5 

0.  831 

0.  906 

6 

0.909 

0.913 

7 

0.909 

0.  941 

8 

0.948 

0.  948 

Table  6.4.  The  Values  of  the  Right  Side  of  (6.14)  (for  A^). 


Since  C.(p.k)  -  1,  the  accuracy  of  the  eigenvalue  approximation  can  be 


read  off  from  Figures  6.3,  6.5,  and  6.6,  changing  the  scale  of  the  error  by 
squaring  -  a  107.  approximation  error  corresponding  to  a  V/.  eigenvalue  error,  a 
17.  approximation  error  corresponding  to  a  .017.  eigenvalue  error,  etc. 


Figure  6.7  shows  the  relative  error  in  A  (p,k)  as  a  function  of  the 


number  of  degrees  of  freedom  N(p,k).  The  polynomial  degree  p  is  also 
shown.  We  see  once  more  the  typical  S-curve. 


si 


,1 


% 


"*4 


44 


RELA‘ 


2°n  Gv2 


- x  k  =  0 

- o  k=  I 

—  +  k=2 


0.001 


2  4  10  20  40  100  200  400  100 

NUMBER  OF  DEGREES  OF  FREEDOM  N 

Figure  6.7.  Relative  Error  in  A  (p,k). 

From  Figure  6.7  we  note  that  the  smallest  eigenvalue  A^  is  approximated  with 
an  accuracy  of  <  3%  for  N  =  2  (corresponding  to  k  =  0  and  p  =  2 ) .  From 
Figures  6.5  and  6.6  and  Table  6.1  we  see  that  A  and  A  are  approximated 


with  this  accuracy  for  N  =  9  (corresponding  to  k  =  0  and  p  =  4).  These 


results  illustrate  a  feature  of  eigenvalue  approximation:  reasonable  accuracy 
for  the  low  eigenvalue  can  be  achieved  relatively  cheaply. 


Accuracy  of  the  Approximate  Eigenfunctions 

Now  we  turn  to  a  discussion  of  the  accuracy  of  the  finite  element  approx¬ 
imation  of  the  eigenfunctions.  As  in  the  discussion  of  the  accuracy  of  the 
eigenvalues  we  define  (cf.  (5.101) 


(6. 15) 


II  u ,-u  ( p,  k )  ||  p/|l  u  || 

D.(p.k)  =  1  *  E - Ll. 

i  K  e. (p,k) 


From  Theorem  5. 1  we  have 
(6. 16) 


1  <  D. (p,k)  <  1  +  0(p 


-4/3, 


Table  6.5  presents  the  values  of  D.(p,k)  for  various  values  of  p  and  k. 


1. 060 


1.062  I  1.062 


1.004  1.015  1.014 


1.001  I  1.000  1.001 


1.014 


1 . 003  1 . 009  1 . 008  1 . 009 

1. 003  1.003  1. 004  1. 004 

1.002  1.001  1.001  1.001 


1. 001 


1.001  |  1.000  1.000  1.000 

1.000  I  1.000  1.000  1.000 


1. 041 
1. 049 
1. 009 
1.006 


1.116  1.115 


1.038 
1.047 
1. 015 
1.001 


1. 038 


1.015 


1.  115 
1. 038 


1.047  1.047 


1. 016 


1.002  1.002 


1.001  1.001  1.002  1.002 

1.000  1.000  1.000  1.000 


1.000  1.000 


1 . 000 


1 . 000 


1. 414 


1.163  1.166  1.166 
1.039  1.040  1.040 


1.041 


1.041  1.041 


1.002  1.019  1.019 


1.002 

1.001 


1.004 

1.004 


1. 019 


1.004  1.004 

1 . 004  1 . 004 


1.001  1.001 


1.001 


1.001 


1.000  1.000  1.000  1.000 


Table  6.5.  The  Values  of  D.(p,k). 

l 

We  clearly  see  that  D^(p,k)  is  very  nearly  1  for  the  entire  range  of  p 
and  k.  Thus  the  accuracy  of  eigenfunction  approximation  can  be  read  off  from 
Figures  6.3,  6.5,  and  6.6.  We  also  see  that  D  >  1,  D.  converges  to  1, 


but  not  monotonical ly,  and  the  convergence  appears  to  be  better  than 


MiSI 


We  see  that  the  values  of  Q^(p,k)  are  nearly  independent  of  p,  which  is 
what  we  would  expect  from  (6.17).  Note  that  the  first  factor  on  the  right 
side  of  (6.14)  is  l-[Q^(p,k)p  1  ]2 . 

Regarding  eigenfunction  error,  (5.12)  implies 

II  u.  -  u.  (p,  k)  ||  b  <  Cp  V.  -ui(p,k)||E, 

To  illustrate  this  result  we  define 


(6.20) 


p||u.— u.  ( p,  k)  || b/|| u.  ||b 
Qi(p,k)  =  ||ui-ui  (p.  k)  ||E/||uillE* 


Table  6.7  presents  values  of  Q.(p,k)  for  k  =  0.  The  results  are  seen  to  be 

similar  to  those  in  Table  6.6.  We  can  thus  conclude  that  the  errors 

u.  -u.  (p.k)  and  u.  -E(p,k)u.  have  very  similar  behavior.  This  is  true 
ii  i  i 

in  either  of  the  norms  11*11.  or  11*11-  and  for  rough  as  well  as  smooth 

D  L 

eigenfunctions,  meshes  with  distorted  as  well  as  undistorted  elements,  and  for 
the  entire  range  of  p  and  k  that  was  considered. 


V 

p  \ 

A1 

*2 

X3 

1 

2 

0.690 

1.  193 

- 

3 

0.  722 

1.934 

2. 567 

4 

0.980 

1. 950 

1.709 

5 

0.964 

1.878 

2.  136 

6 

0.953 

1.382 

1.850 

7 

0.  841 

1.304 

2.  158 

Table  6.7.  Values  of  Q,(p,0). 


Table  6.8  shows  the  values  of 


where  ||u.||,  =  ||u.(p,k)||,  =  1,  for  k  =  0  and  1  and  i  =  1.  We  see  (as  we 
lb  l  b 

should  expect)  that  the  first  of  these  two  values  is  the  slightly  smaller  than 
the  second. 


p 

l|Ui  — Ui  (  P ,  0  )  II  E 

II  ut— Ui  (  P,  0  )  II  E 

llui-u^p,  1 )  ||E 

llui— Ui  (p,  1 )  ||E 

II  Ui  II  E 

II  Ui  II  E 

II Ui  ||  e 

II St  Be 

1 

- 

- 

0.  666 

0.  800 

2 

0.  174 

0.  176 

0.339 

0.  352 

3 

0.  149 

0.  150 

0.235 

0.  239 

4 

0.  132 

0.  133 

0.088 

0.089 

5 

0.  109 

0.  109 

0.  032 

0.033 

6 

0.  086 

0.086 

0.  026 

0.026 

7 

0.071 

0.072 

0.020 

0.  021 

8 

0.062 

0.  062 

0.  017 

0.  017 

Table  6.8.  Comparison  of  the  Errors  in  u  (p,k)  and 
u^p.k)  in  ||  *  ||  E- 


Table  6.9  presents  the  corresponding  errors  in 


1  b’ 


l .  e.  , 


I  ui_u1  (P>  k )  II b  || u1~u1  (p,  k)  || b 

-  and  - 


u 


1  b 


T'b 


for  k  =  0  and  1.  The  table  shows  that  the  second  of  these  quantities  is 


p 

Huj-ut  (p,0)||b 

C  1 

w-» 

1 

C  1 

V 

o 

O' 

llui-ujlp,  1  )  Hb 

|| u1-u1  ( p,  1  )  ||b 

luillb 

llUiUb 

llutllb 

II  U,  Hb 

1 

- 

- 

0. 3175 

0.2831 

2 

0.0601 

0.0590 

0. 1252 

0. 1175 

3 

0.0359 

0.0345 

0. 0801 

0.0771 

4 

0.0323 

0.0314 

0.0029 

0.0226 

5 

0.0210 

0.0202 

0.0035 

0.0034 

6 

0.0136 

0.0131 

0.0020 

0.0020 

7 

0.0079 

0. 0075 

0.0005 

0.0004 

8 

0. 0068 

0. 0064 

0.0004 

0. 0004 

Table  6.9.  Comparison  of  the  Errors  u^lp.k)  and 


u  (p,k)  in 


In  general,  to  measure  the  error  between  exact  and  approximate  eigen¬ 
functions  we  can  associate  the  exact  and  approximate  eigenfunctions  together 


in  various  ways.  For  example,  we  can  associate  u^.  and  u.(p,k)  or  u.  ar 


u  (p,k).  A  third  possibility  is  to  associate  u^lp.k)  with  u^  =  /3u.,  wher 


13  is  chosen  so  that  ||u.  -u.(p,k)||_  is  minimal.  This  choice  is  of  interest 

l  l  E 


since  it  also  minimizes 


I /3ui -ui  ( p,  k )  ||  b 


II Pu.  -u.  (p,  k )  I! E 


Table  6. 10  shows  the  values  of 


I  u1-u1  ( p,  k )  II E  II Uj-Uj  ( p,  k )  II  b 

- and  - . 


Note  that  II  u  1 1| ^  <  Hu^. 


riv.v.v.VwWA'.VA'  v  hjih.*-  o  w.v.v.v.v.v'. 


VI  V*  T*»  ' 


^77777777  v-¥v'.v/j  ',-w  '/  *7’  .'.v.v.v.v' 


p 

lluj-Ui  (p,  0)  ||E 

HUj-U!  (p,0)||b 

llUi-u^p,  1 )  || E 

HUj-Uj  (p,  1  )||b 

l|u||  E 

|ut|b 

II  U1  II  E 

llu||  b 

1 

- 

- 

0.  808 

0. 2912 

2 

0.  176 

0.0590 

0.354 

0. 1181 

3 

0.  150 

0.0345 

0.  240 

0. 0773 

4 

0.  133 

0.0314 

0.089 

0.0226 

5 

0.  109 

0.0202 

0.032 

0. 0034 

6 

0.086 

0.0131 

0.026 

0.0020 

7 

0.071 

0.0075 

0.  021 

0.0004  ! 

8 

0.062 

0.0064 

0.017 

0.0004 

Table  6.10.  Comparison  of  the  Errors  u^-u  (p,k) 


in 


and 


'E  “““  "  'b‘ 

We  see  that  these  different  associations  of  exact  and  approximate  eigenfunc¬ 
tions  lead  to  very  similar  values  for  the  errors. 


A  Posteriori  Estimates 


Let  (A,u)  =  (A^.u^  be  a  fixed  eigenpair  of  (6.1)  and  let  (A(p).u(p)) 


=  (A. (p, k) , u^ (p, k) )  be  its  approximation  by  the  finite  element  method.  We 


suppress  both  i  and  k  in  this  notation.  Denote  by  z^P^  €  the 


exact  solution  of 


(6.2i: 


-Az 


=  A(p)u(p)  in  Q 


l  pi  n  _ 

z  =0  on  f, 


and  then  denote  by  z 


[p] 


is 


the  finite  element  approximation  of  z 
P) 


•  00.  Let 


the  finite  element  solution  of  (6.21)  using  elements 

u(p) 

H 

\  rA  1  f  r>  1 

and 

Hz' 


of  degree  q  >  p  on  the  same  mesh.  Obviously  z^P  =  u(p) 


Ipl 


|zlpJ  -zlpJ 


z^  -  u(p)  II  and  e 

t,  CJ 


[Pi  =  njpl 


q 

u  (  p )  ii 


■  0  as 


It  is  shown  in  [17]  that 


52 


[6.22: 


|A(p)~A|  <  (c[pl )2( l+o(l) ). 
A 


so  that  if  we  have  a  reliable  a  posteriori  estimate  for  c 


we  then  have  a 


reliable  a  posteriori  estimate  for  -A— — — — .  There  are  a  number  of  ways  to 


A 

,[p) 


obtain  an  a  posteriori  estimate  for  elf".  For  a  survey  of  relevant 
literature  we  refer  to  [18].  We  now  obtain  such  an  estimate.  Note  that 


(6.23: 


clqp]  =  (Hz^1!!2  -  ||u(p)||2)1/2. 


and  for  q  >  p, 


f  [pH 2  ..  [p]  ,  ,  1, 2  „  Ip]  Ip]  .,2 

\cq  =  II z  -u(p)||E-||z  r  --Z.J  || r. 


'q  11 E' 


.tpl 


We  thus  easily  see  that  cq  is  increasing  with  q  and  that 


(6.24) 


1  im 
q->oo 


.IP' 


.[Pi 


=  llz^1  -  u(p)  I 


retp; 


Hence 

(6.25) 


c[p[a,£.[p[>  for  q  large; 

q 


this  is  the  desired  a  posteriori  estimate  for  c 

2 


P) 


(6.26) 


|  A(p)-A| 


■  Kpl)  ■ 


From  (6.22)  we  have 

for  1  «  p  «  q. 


To  check  the  quality  of  the  approximation  in  (6.26)  we  consider  the  quantity 


(6.27) 


,[p]  _  ( A ( p ) -A ) /A 


Kpl)' 


Jp] 


In  Tables  6.11  and  6.12  the  values  of  are  presented  for  various  values 


of  p  and  q  and  for  0-  and  1-layer  meshes  for  the  first  eigenvalue.  Q 
is  obtained  by  an  extrapolation  procedure. 


[Pi 
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